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In isotropic helical turbulence, a new single helical model is suggested for large eddy simulation. Based on 
the Kolmogrov's hypotheses, the helical model is proposed according to the balance of helicity dissipation 
and the average of helicity flux across the inertial range, and the helical model is a kind of eddy viscosity 
model. The coefficient of the helical model is constant and can be also determined dynamically. Numerical 
simulations of forced and decaying isotropic helical turbulence demonstrate that the helical model predicts 
the energy and helicity evolution well. The statistical character of the helical model is closer to the DNS 
results. In a priori test, the energy and helicity dissipations of the helical model also show a good results. In 
general, the helical model has the advantage contrast with the dynamic Smogorinsky and mixed models. 



I. INTRODUCTION 

In turbulence flow, helicity is an important physical 
quantity. It is widespread in the motions of the atmo- 
sphere, ocean circulation and other natural phenomena, 
and also found in leading edge and trailing vortices shed 
from wings and slender bodies^—. Helicity can be de- 
fined as, h = u ■ u> , where u and oj are the velocity and 
vorticity of the turbulence flow respectively, and h is a 
pscudoscalar quantity. Similar with the status of energy 
in the dynamics of ideal fluids, helicity has the character 
of inviscid invariance. This physical property determines 
that helicity is an important quantity in turbulence re- 
search. 

Recently, the researches on helical turbulence have de- 
veloped greatly in theoretical research, experiment and 
numerical simulation. Based on the helical decomposi- 
tion of velocity, the mechanism of existing a joint for- 
ward cascade of energy and helicity has been explained 
in theory^. Cascades existing in helical turbulence have 
space scale and time scaled— , and the researches showed 
that the existing space scale of helicity cascade was large 
than energy cascade. In the inertial range, the joint cas- 
cade of energy and helicity was dominated by the energy 
cascade time scale in the low wave number and by the 
helicity cascade time scale in high wave number. Using 
direct numerical simulation of isotropic turbulence, en- 
ergy and helicity flux were studied, the research showed 
that helicity flux is more intermittent than the energy 
flux and the spatial structure were much finer—. 

Large eddy simulation (LES), as a important method, 
has been used to research helical turbulence. Several 
kinds of SGS models have been proposed so far, such 
as eddy-viscosity model, dynamic model, vortex model, 
subgrid-scalc estimation model, and even the developing 
constrained SGS modeU^ - — . Although so many types 
of SGS model, there are few models on account of the 
character of helical turbulence. Y. Li et al. proposed a 
two-term dynamic mixed models (DSH) to do large eddy 
simulations of helical turbulence^. Compared the LES 
results of the DSH model with other traditional mod- 
els, the improvement of DSH model was not so reamark- 
able.In this paper, we first use the balance of helicity 



dissipation and helicity flux across the inertial range to 
get the eddy viscosity, then deduce a new helical model 
(SR). Using SR model and dynamic SR model (DSR)to 
do large eddy simulation of isotropic helical turbulence, 
we find that SR and DSR models can predict energy and 
helicity well. Contrast with dynamic Smogorinsky model 
(DSM) and DSH model, the new helical model is a simple 
model and its results of LES has get an obvious improve- 
ment. 



II. THEORETICAL ANALYSIS AND SGS MODEL 
CONSTRUCTION 

In isotropic helical turbulence, the forced helicity con- 
trol equation can be get from the N-S equation, as below 

8 t h + dj ( Uj h) = djflj - AvSijlhj + 2fiUi, (1) 

where Sy = | (djUi + diUj) is the strain rate tensor and 
Rij = -i (djUJi + diLUj) is the symmetric vorticity gradient 
tensor. f { is the force, and flj is the flux term, which can 
be expressed as 

P 1 

% = h -UiUiUjj + 2v(uiRij + uJiSij) - SjkiUkfi- (2) 

9 * 

For LES, the resolved helicity can be defined as h A = 
u u), and we can get the control function of h A — through 
filtering the Eq.(l) at scale A as 

d t h A + dj (ujh A ) = djQj - U H - \i'S,,l>',, + 2ffii, (3) 



where Tin 



-2rijRij is the SGS helicity dissipation 



rate, and Ty- = UiUj — UiUj is the SGS stress. 

The equation for the resolved helicity in statistically 
stationary forced isotropic helical turbulence can be ob- 
tained through taking the ensemble average of Eq.(3) as 



2(f i u l ) = -2(ryi? jj -)+4^(5yi? i3 -). 



(4) 



where 2{fi0Ji) = rj (rj is the total helicity dissipation) 
—2{njRij) is the SGS helicity dissipation, and Au{SijRij) 
is the viscous helicity dissipation at scale A. 
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In the same way, we can also get the equation for the 
resolved energy in statistically stationary forced isotropic 
helical turbulence^from the filtered NS equation as 



(/*«») = -(njSij) + 2u(SijSij}. 



(5) 



where {fiUi} = e (e is the total helicity dissipation), 
— (TjjSij) is the SGS energy dissipation and 2v(SijSij) 
is the viscous energy dissipation at scale A. 

From Eq.(4) and Eq.(5), we can see that the helicity 
dissipation has the similar form and composition as the 
energy dissipation. 

In the incrtial range of isotropic helical turbulence, the 
viscosity of the flow may be ignored, therefor we can get 
the helicity and energy relations^ 



and 



i] = -2{njRij), 



-\ T ijSij). 



(6) 



(7) 



At the same time, the average of the energy flux across 
the inertial range is invariable and equals to the SGS 
energy dissipation, 



(He) = -(njSij), 



(8) 



where IIe is the energy flux across the inertial range. 
Similar with energy, helicity has the same character in 
incrtial subrange, 



(U H ) = ^(TijRij), 



(9) 



where Hh is the helicity flux across the inertial range. 

In large eddy simulation of isotropic helical turbulence, 
we take the eddy- viscosity model as 



■^SijTkk = -2l> t S., 



%3 ! 



(10) 



where vt is the eddy viscosity. 

Based on the assumption of Eq.(9) and Eq.(10), we 
can obtain 

(U H ) = -2(T™ od R tl ) = 2 • (2v T S l3 R l3 ), (11) 

where r, r ™ od is the SGS stress model. 

In this deduction, we consider the eddy viscosity vt as 
a global averaged valued for the time being, thus 

<n ff > = 2/< ; 2.S-,,//,, . (12) 

In isotropic helical turbulence, we have 



and 



{2S tj S l3 ) = / 2k 2 E(k)dk, 
Jo 



2(2S ij R ij ) = / 2k 2 H(k)dk, 



(13) 



(14) 



where E(k) and H(k) are the energy and helicity spectra 
functions, and k c is the cut off wavenumber. 

The energy and helicity spectra E(k) and H(k) have 
the inequalit y 23 ' 24 as 



\H(k)\ < 2kE{k). 



(15) 



In the inertial range, the energy and helicity spectr a^ 
can be choosed as 



and 



E(k) = C K e 2 / 3 k- 5 / 3 , 



H{k) = C H m~ 1/3 k- 5/3 , 



(16) 



(17) 



Ck is Kolmogrov constant and Ch is the coefficient of 
the helicity spectrum. 

Substituting Eq.(17) into Eq.(14), we can get 



From Eq.(15), we can see clearly that 
\H(k)\ 



E{k) 



ak. 



(18) 



(19) 



where a is a dimensionless parameter, its value changes 
with k, and < a < 2. 

When the cut wavenumber is in the inertial range, from 
Eq.(16), Eq.(17) and Eq.(19), we obtain 

C H 



C k 



t 1 k~ 1 \i 1 \ 



(20) 



Taking the absolute value of both side of Eq.(18), and 
putting it into Eq.(20), we can get the expression 

2\(2S ij R ij )\ = f (aCWC|) 1/3 M 2/3 ** /3 . (21) 

Then, from Eq.(6), (9), (11) and (21), we can obtain the 
expression of vt-, 



(22) 



where A = Tr/k c is the filter scale, S r =\(2SijRij)\ 1 / 2 and 
the model coefficient C r = (l) 3/2 n- 5 / 2 (aC K C 2 H y 1 ' 2 . 
Thus, we obtain a new helical SGS model as 

In the LES model, a is set to 2.0 and Ck = 1-6 and 
Ch is confirmed by the DNS data about 1.35. 

In this letter, the eddy viscosity vt of the model may 
be given a local value, and the expression of vt is 



v T = lc r A 5 / 2 S SR , 



(24) 



where Ssr = \2SijRi 3 -] 1 / 2 . Thus we get the new helical 
model discussed in this letter (SR model) as 

r™ od = -Cr^SsRSij. (25) 

where C r « 0.036. The coefficient of helical model C r 
can also be decided dynamically, and then we get the 
DSR model in this letter. 
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III. THE NUMERICAL RESULTS AND ANALYSIS 

In this part, we will give a priori and a posteriori test 
of the LES model, and do some comparison and analysis. 

In order to validate our model, a DNS of three- 
dimensional incompressible homogeneous isotropic tur- 
bulence is introduced here. It solves the forced N-S equa- 
tions using a pseudo spectral code in a cubic box with 
periodic boundary conditions, and the numerical resolu- 
tion is 512 3 . A Guassian random field is the initial flow 
condition, and it has an energy spectrum as 



E (k) = Ak 2 U^k t 



(26) 



where ko = 4.5786 and Uq = 0.715. The whole system is 
maintained by a constant energy input rate e = 0.1 and 
r\ = 0.3 in the first two wave number shells. 

Have constructed the new helical models SR and DSR, 
we first test the validity of the models a priori. In the 
inertial range, we have assumed the invariance of the en- 
ergy and helicity SGS dissipations. In such a precondi- 
tion, the invariance of f(S) and h(S) must hold, 



/(*) = (S 



2.") 



SsiiSijSij), 



(27) 



and 



h(S) = (S^SsRSijRij), (28) 
where 5 is the length scale varying in the inertial range. 
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FIG. 1. (a)/(<5),(b)/i(<5) distribute with S/( for a priori. ( is 
the Kolmogrov length scale. 

In Fig.l, we show f(S) and h(S) as a function of S/(, 
and £ is the Kolmogorov scale. It is easy to see that the 



values of the two functions in inertial range are almost 
constant, and it demonstrates the scale-invariance and 
the reasonableness of the model. 

We choose four LES models to do some analysis and 
comparison here, the SR model, the DSR model, the 
dynamic Smogorinsky model (DSM) and the dynamic 
mixed helical model (DSH)2I. 




FIG. 2. Energy spectra for steady isotropic turbulence. Solid 
line: DNS; dashed line: SR; line with square: DSR; dash-dot- 
dotted line: DSM; line with delta: DSH. 
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FIG. 3. Helicity spectra for steady isotropic turbulence. Solid 
line: DNS; dashed line: SR; line with square: DSR; dash-dot- 
dotted line: DSM; line with delta: DSH. 

Fig. 2 and Fig. 3 display the energy and helicity spectra 
of DNS and different models for steady isotopic turbu- 
lence respectively. The bolt solid line is for DNS spectra, 
and the other lines are for the four models. It is easy to 
see that the SR and DSR show similar trends and pre- 
dict the energy and helicity spectra better than DSM and 
DSH. In Fig.2, DSM and DSH underestimate the energy 
spectra near the cutoff wave-number and overestimate 
at the middle range of the wave-number. In Fig. 3, DSM 
and DSH overestimate the helicity spectra at the middle 
range of the wave-number, and DSH underestimates it 
seriously near the cutoff wave-number. 
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FIG. 4. Energy spectra for decaying isotropic turbulence 
(aposteriori), at t = 0, 6ro, and 12ro , where to is the ini- 
tial large eddy turnover time scale. Solid line: DNS; dashed 
line: (a) SR, (b) DSR, (c) DSM, (d) DSH . 



FIG. 5. Helicity spectra for decaying isotropic turbulence 
(aposteriori), at t = 0,6to, and 12tq , where to is the initial 
large eddy turnover time scale.. Solid line: DNS; dashed line: 
(a) SR, (b) DSR, (c) DSM, (d) DSH . 



In Fig. 4 and Fig. 5 wc show the time evolution of the 
DNS and four models' energy and helicity spectra for a 
decay problem starting from a fully developed statistical 
steady state respectively, where t = 0, 6t and 12r . In 
Fig. 4, similar with the steady prbolem, the SR and DSR 
also show similar trends and predict the energy spectra 
very well. DSM and DSH overestimate the energy spectra 
at the middle range of the wave-number and underesti- 



mate it near the cutoff wave-number. Fig. 5 show us the 
graphics different from Fig. 3 greatly. The helicity spec- 
tra have large fluctuation, because it is a free decaying 
course and helicity is a pseudoscalar quantity , the two 
factors cause such phenomenon. In Fig. 5 we can still sec 
SR and DSR pridict the helicity spectra a little better 
than the other two models. 

Fig. 6 and Fig. 7 show us the SGS energy and helicity 
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FIG. 6. The SGS energy dissipation (steady) distribute with 
5/C for a priori. Bold solid line: DNS; line with Square: SR; 
line with delta: DSR; line with diamond: DSM; line with 
circle: DSH. 




FIG. 7. The SGS helicity dissipation (steady) distribute with 
5/C for a priori. Bold solid line: DNS; line with Square: SR; 
line with delta: DSR; line with diamond: DSM; line with 
circle: DSH. 



FIG. 8. The decay of SGS energy dissipation with t/ro (to 
is the initial large eddy turnover time scale) from a fully de- 
veloped steady state for a priori at the filter scale A . Bold 
solid line: DNS; line with Square: SR; line with delta: DSR; 
line with diamond: DSM; line with circle: DSH. 




FIG. 9. The decay of SGS helicity dissipation with t/ro from 
a fully developed steady state for a priori at the filter scale A 
. Bold solid line: DNS; line with Square: SR; line with delta: 
DSR; line with diamond: DSM; line with circle: DSH. 



dissipations of DNS and the four models distribute with 
8 / Q for a priori respectively, and they reflect case of the 
full developed steady turbulence. It is obvious from in 
the inertial range that the SGS energy and helicity dissi- 
pations from SR are closer to those of the DNS than other 
models, and specially have the similar change trends with 
DNS. 

In Fig. 8 and Fig. 9, we show the decay of SGS energy 
and helicity dissipation with £/to(to is the initial large 
eddy turnover time scale) from a fully developed steady 
state for a priori at the filter scale A respectively. It is 
clear in Fig. 8 that the SGS energy dissipations from SR 
and DSR are closer to the DNS result than DSM and 
DSH, and also have similar change trend. And also we 
can see in Fig. 9 that the SGS helicity dissipations from 
SR, DSR and DSM are closer to the DNS result than 
DSH. 

In Fig. 10, we display the distribution of high-order ve- 



locity increment with A, where r is the separation dis- 
tance and A is the filter scale. Fig. 10 (a), (b) and (c) de- 
note the fourth-order, sixth-order and eighth-order struc- 
ture functions of the longitudinal velocity increment of 
DNS and the four models distribute with separation dis- 
tance. We can from Fig. 10 that the result from SR is clos- 
est to the DNS result very well in the range of r/A > 3, 
and the behavior of DSR is a little better than DSH and 
DSM. While in the range of r/A < 3, anyone of the 
models have still not given a rather good result. 



IV. CONCLUSIONS 

In this letter, we construct a new helical SGS stress 
model for large eddy simulation of the isotropic helical 
turbulence. Different other eddy viscosity model, we con- 
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FIG. 10. The distribution of high-order moments of longitu- 
dinal velocity increment with the separation distance r, where 
A is the filter scale, (a) Sf , (b) Se , (c) S£ 



duct the helical model based on the balance of helicity 
dissipation and the average of helicity flux across the in- 
ertial range in helical turbulence. Then we have tested 
the scale-invariance of the model for a priori in inertial 
range, and make sure the validity of the helical model. 
We have given a constant coefficient helical model and a 
dynamic helical model here. 

Have tested from a priori and a postoriori, the helical 
model is confirmed to predict energy and helicity spectra 
accurately, and also gives rather good simulating results 
in energy and helicity dissipation, et aLBesides the ad- 
vantages above all, the helical model is a single model 
and can improve the computing efficiency greatly. As a 
SGS model on account of helical turbulence, it can also 
apply to the rotational turbulence. 
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